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Abstract 

Hirota's bilinear method ("direct method") has been very effective in construct- 
ing sohton solutions to many integrable equations. The construction of one- and 
two-soliton solutions is possible even for non-integrable bilinear equations, but the 
existence of a generic three-soliton solution imposes severe constraints and is in 
fact equivalent to integr ability. This property has been used before in searching 
for integrable partial differential equations (PDE), and in this paper we apply it to 
two dimensional partial difference equations (PAE) defined on a 3 x 3 stencil. We 
also discuss how the obtained equations are related to projections and limits of the 
three-dimensional master equations of Hirota and Miwa, and find that sometimes a 
singular limit is needed. 



1 Introduction 



Integrable dynamics is fundamentally associated with regularity and predictability, but 
if the equation is nonlinear this does not imply that the behavior is simple, quite the 
contrary. 

There are in fact many interesting mathematical properties that are associated with 
nonlinear integrable PDE's (soliton equations), such as the existence of sufficient number 
of conserved quantities and symmetries, the ability to write the PDE as the consistency 
condition of two linear equations (the Lax pair, with spectral parameter), the existence 
of multi-soliton solutions, etc.[l]. Any of these characterizing properties can also be used 
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to search for integrable PDE's and it is remarkable that in the end each criterion gives 
essentially the same set of equations. 

In recent years the study of integrable partial difference equations (PAE) has pro- 
gressed rapidly [2]. The above integrability indicators have also been applied to PAE's, 
but in the discrete case there are also other criteria, specific to discrete equations, such 
as singularity confinement [3], the so called "Consistency- Around-the-Cube" property [4], 
and slow growth of complexity of the iterates (algebraic entropy) [5]. 

In this paper we consider the existence of three-soliton solutions for PAE's as an indica- 
tor of integrability We use Hirota's bilinear method, which is well suited for constructing 
soliton solutions. In Section 2 we first describe how the method works in searching for 
integrable PDE's and then formulate the search problem for PAE's. In Section 3 we 
give the results for a particular class of two-dimensional lattice equations, namely those 
defined on a 3 x 3 stencil. It is remarkable that all of the integrable equations obtained 
from this class can in fact be described as projections and limits of the well known three- 
dimensional bilinear master-equations or Hirota [6] or Miwa [7], although sometimes a 
singular limit is needed. This is shown in Section 4. 

2 Hirota's direct method 
2.1 The continuum case 

In 1971 R. Hirota proposed a direct method for constructing soliton solutions [8]. A key 
ingredient in this method is a transformation to new dependent variables so that in these 
new variables a soliton solution is a finite sum of exponentials. 

2.1.1 Definition 

For the prototypical example of KdV equation 

ut + u^cxx + 6wwx = 0, 

one applies the transformation 

u = 29^ log F, 

and after integrating the resulting equation once one gets a bilinear equation 

{Dt + D,Dt)F-F^O, (1) 

were the Hirota derivative operator D is defined by 

D2f-9 = {d^r-d.XfixMx,)l^^^^^^ (2) 
= d^f{x + y)g{x-y)\^^^. 

Thus D acts on a product similarly to the Leibniz rule, except for a crucial change of sign. 
In particular f-g — {—l)'^D"g-f. Equations in which derivatives appear only through 
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the above D-operaior are said to be in Hirota bilinear form. (For further properties of 
the D operator, see the book [9] and Appendix I of [10].) 

The first step in constructing sohtons solutions is therefore to find the dependent vari- 
able transformation that leads to a Hirota bilinear form. Unfortunately the bilinearization 
of a given nonlinear PDE is not easy. One cannot know before hand how many dependent 
or even independent variables one need to properly bilinearize a given equation.^ 

2.1.2 Gauge invar iance 

At this point we would like to emphasize that equations written in Hirota bilinear form can 
also be characterized by gauge invariance. Consider the gauge transformation {f,g,...) i->- 
{e^'^f, e^'^g, . . .), then the condition for gauge invariance is 

P{D){e^-^f) . {e^-^g) = e^'-^P{D)f ■ g. (3) 

(Here we use vector notation: D = {D.^, Dy, D^, . . .), x = {x,y, z, ■ ■ ■), p = {Px,Py,Pz, ■ ■ ■ ) 
with inner product x ■ p = Px x + Py y + ■ ■ ■ .) Equations written in Hirota bilinear form 
have this property, and conversely, any quadratic equation having gauge invariance of the 
above type can be written in Hirota bilinear form [11]. 

Gauge invariance allows generalizations of the Hirota bilinear form to other circum- 
stances, such as higher multi- linearity [11] and discrete equations. 

2.1.3 ISS and 2SS 

The existence of a bilinear form does not guarantee the existence of A^-soliton solutions, 
but one- and two-soliton solutions are often easy to construct. As an example take the 
bilinear equation 

P{D)F • F = (4) 

where P is some even polynomial in D (such as (1) in the KdV case). 

In order to start the expansion we need a vacuum or background solution, which in 
this case is F = 1 (which in the KdV case corresponds to u = 0). This works, provided 

that P(0) = 0. 

Starting with the vacuum F = 1 one builds a one-soliton solution (ISS) "perturba- 
tively" using exponentials of hnear functions, which are called plane wave factors (PWF): 

F — 1 + ee^\ where rji — x ■ Pi + rj^, with constant 77°. (5) 

Here e is the formal expansion parameter, and upon substituting this to (4) one finds the 
constraint that the parameters Pi must satisfy the dispersion relation (DR) 

P{Pi) = 0. (6) 

^In the literature one sometimes finds examples of over-bilinearization, in which a multi-linear equation 
is broken down into too many bilinear equations. One case study can be found in [12]. 
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Continuing in this manner one finds that the generic equation (4) also has two-sohton 
solutions of the form 

F = 1 + e + e e''^ + ^126"^+"^ , (7) 

where 

Here each pi is restricted by the DR (6), that already appeared in constructing the ISS, 
but there are no additional restrictions on them. Aij is called the phase factor. Note that 
Aij = Aji due to P being even. 

In the above constructions crucial role is played by the minus sign in the definition 
(2) of the Hirota derivative, which implies 

P{D)e^-P ■ e^-P' = P{p- p')e^-^+P'\ (9) 

In practice this means that the first and the last terms of expressions like the ISS (5) and 
2SS (7) automatically satisfy the equation (if P(0) = 0). 

The above construction shows a level of partial integrability: we can have elastic scat- 
tering of two solitons, for any dispersion relation, if the nonlinearity is suitable (namely 
if it arises from a bilinear equation as above) . 



2.1.4 The 3SS-condition 

Definition 1 (Hirota integrability) 

Assume that the bilinear equation (4) has ISS's of the form 

F^l + ee^^, r]j^x-pj + r]^, 

where the parameters pj are only restricted by the DR P{pj) = 0. The equation is said to 
be Hirota integrable if it has NSS's of the form 

N 

F — 1 + e e''J + (finite number of h.o. terms in e) 

without any further conditions on the parameters pj of the individual solitons. 

The assumption that there are no further restrictions is important, because any equation 
has multi-soliton solutions for some restricted set of parameters. 

If we now apply this principle to the 3SS we start with the ansatz 

F = 1 + e e''! + e e''^ + e e''^ 

Aue'^'+^' + ^236"^+"^ + Agie^^+^i + Al2A23Al3e''^+''^+''^ (10) 

which is fixed by the requirement that if in a NSS any soliton goes far away, the rest should 
look hke the (N-l)SS. (In practice "going away" means either 6*"= — >■ or e''* — >■ 00.) This 
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means that there is no freedom left: parameters are restricted only by the dispersion 
relation (6) and the phase factors Aij were given already in (8). Thus the existence of 
three-soliton solutions is not automatic, but rather it imposes severe requirements on the 
polynomial P. 

The above definition of integrability can be used as a method for searching for new 
integrable equations. In practice it means that one should find polynomials P that verify 
the "three-soliton-condition" (3SC), which follows when the ansatz (10) is substituted into 
(4), namely 



P{criPl + (72P2 + (73P3)P((TiPi - (T2P2) 

XP((72P2 - Cr^P^)P{(TlPl - CTsPs) = 0, 



=±1 



(11) 



on the manifold P{pi) — 0, Vi //the 3SC is satisfied then the equation most hkely has 
NSS of the form 



J2 exp 

we{o,i} 



(N) N 
i>j 1=1 



where exp(aij) = Aij. 



(12) 



2.1.5 Continuous searches 

Condition (11) can be easily converted for symbolic computation, for example in REDUCE 
[13] one can use the LET-rule for the dispersion relations (which in REDUCE is evaluated 
to the end). In 1987-1988 one of the present authors (JH) used this method to search for 
new integrable equations of KdV, mKdV, sG and nlS type [14, 15], for which conditions 
similar to (11) can be derived. Among the (very few) new equations found in this way 
were 

{DlDt + aDl + D,Dy)F • F = 0, (13) 
{Dl - D^Dl + aDl + hD^Dt + cDl)F • F = 0, (14) 

and 

/ {iaDl + 3D,Dy - 2tDt + c) G ■ F = 0, 

\ [a{a^Di-3Dl + 4aD,Dt) + bDl]F-F = \G\\ ^^^^ 

The pair of equations (15) is of the nonlinear Shrodingcr type and in fact combines the 
two most important (2 + 1) -dimensional equations, Davey-Stewartson and Kadomtsev- 
Petviashvili equations, but only their DSII and KPI variants. 



2.2 Hirota's bilinear formalism for lattice equations 

We now turn to discrete lattice equations and as a model we use the continuous situation 
described above. The various differences between continuous and discrete models (such 
as the lack of Leibniz rule in the latter) make things now much more difficult. We will 
only consider PAE's defined on the 2-dimensional cartesian lattice. 
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2.2.1 Discrete Hirota bilinecir form 

Before we can proceed to equations we must have a definition of the discrete Hirota 
bihnear form. The guiding property is still the gauge-invariance (3), and thus assuming 
we have lattice functions m), we require gauge-invariance under 

/,(n,m) ^ /;(n,m) = A"5™/,(n,m). (16) 

Guided by the requirement of invariance under (16) we have 

Definition 2 (Discrete Hirota bilinear form) If the equation can be written as 

Yl {n + u^,m + //+) Qj {n + iy-,m + nj) ^0 
j 

where the index sums /I'^' + iiJ — /i^, v'^ + vj — do not depend on j, we say it has 
discrete Hirota bilinear form. 

Note that 

e'^^^ fix) ■ g{x) = e''^^^-^^') f{x)g{x%,^^ = f{x + a)g{x - a) . 

Therefore the discrete version of P{D)f ■ f — should be an even function built up from 
exponentials of the Hirota derivative. For example a sum containing terms like cosh(£)^) 
or sinh(L'a;) sinh(L't) etc., is acceptable. 



2.2.2 The Hirota and Miwa equations 

The first results on discrete bilinear soliton equations were obtained by Hirota in a series 
of papers in 1977 [16]. Hirota's aim was to discretize known soliton equations so that their 
multi-soliton solutions would still have more or less the same form. In this way he was able 
to derive integrable 2D lattice equations for KdV, Toda, and sine-Gordon equations. The 
culmination of this series was the "Discrete Analogue of a Generalized Toda Equation" 
(DAGTE) or "Hirota-equation" [6], which in the fully 3D version (discrete bilinear KP 
equation) can be written as [7] 

a{b — c) Tn+l,m,kTn,m+l,k+l + b{c — a) Tn,m+l,kTn+l,m,k+l + c(a — b) Tn,m,k+lTn+l,m+l,k = 0. 

(17a) 

This equation has soliton solutions of the form (12), but with discretized PWF's given by 
[7, 17] 



i^Vfi^VV^^l , A^ = ^ii^iP4^^i^. (17b) 

1-apiJ \l-bpij \l-cpij [Pi - Qj){Qi - Pj) 

In the early 1980's the continuum bilinear approach was developed to a complete 
theory by Kyoto school, and it was very soon applied to the discrete case as well [7, 17]. 
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We note in particular the generalization of Hirota's equation (17a) by Miwa (the discrete 
BKP equation) [7] 

(a + b){a + C){b - c) rn+l,rn,k'rn,m+l,k+l + {b + c){b + a){c - a) Tn,m+l,krn+l,m,k+l 
+ (c + a) (c + b) (a - b) Tn,m,k+irn+l,m+l,k + (a - 6) (6 - c) (c - a) Tn+l,m+l,k+lTn,m,k = 0. 

(18a) 

Also this equation has soliton solutions of the form (12), except that now [7, 17] 

^ / (I - ap,){l - aq,) y / (I - bp,){l - 6g,) \- / (I - cp.)(l - cq,) y 

\{l + ap,){l + aq,)J \(1 + bp,){l + bq^) J \{1 + cpi){l + cq^) J ' 

(Pi + Pj) {Pi + Qj) {Qi + Pj){Qi + Qj) ' 
2.2.3 Application of Hirota's method to some PAE's 

Hirota's direct method has also been applied to construction of soliton solutions to non- 
linear lattice equations. There the first problem is to construct the dependent variable 
transform that takes one from the nonlinear form to the Hirota bilinear form with gauge 
invariance. This approach was used for most of the equations in the ABS list in [18, 19]. 
For example, in the case of the HI equation 

m = {u-u){u-u) + {a^ -b'^) = 0, (19) 

one is led (after studying the ISS and 2SS) to the dependent variable transformation 

ulT^^an + bm + X-'-p^, (20) 



fn,Tn 

which yields the bilinearization 

HI = - + (a - b)f}] + {a + b)ff\/{ffff) + {a' - b'), 

where 

:Ki = gf-gf + {a-b){ff-f7) = 0, (21a) 
5{2 = 9f-df + (a + b){ff-ff)^0. (21b) 

3 Searching for integrable bilinear lattice equations 

Following the continuous example we would like to search for integrable lattice equations 
by the three-soliton condition. Since we will here only consider 1-component equations, 
we can still take (4) as the basic class of equations, but now P would have to be a sum 
of exponentials. 
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Figure 1: Possible stencils for 2D bilinear equations. Case a) only contains 1 equation 
(22), Case b) is analyzed in this paper, Case d) and similar models are analyzed in [20]. 

3.1 Organizational questions 

The first question is: How to organize such a search? In the continuous case we used the 
total degree in the bilinear derivative and then the degrees in various coordinates. The 
problem in the discrete case is that there are so many ways to discretize a derivative and 
thus there will be many more cases. 

One way to proceed is by the radius of points needed, see Figure 1. For sub-figure la) 
the corresponding bilinear equation will be 

Cl fn,mfn+l,m+l ~l~ C2 fn+l,mfn,m+l — 

0, (22) 

since these are the only combinations of fn+u,m+^J.fn+u',m+^J.' with non-negative n, /x', u, v' 
and sum n -\- ^' — l,v -\- u' — 1. Since we want / = 1 to be a solution we must have 
Ci + C2 = and then there are no parameters, except the trivial scaling. 

In this paper we will consider the arrangement given in sub-figure lb) which corre- 
sponds to the difference equation 

5*8 fn+l,m+l fn—l,m—l + 96 fn+ l,m—l fn—l,Tn+l 
+9'7 /n+l,m /n-l,m + ^6 fn ,m+l /n,m— 1 

+ ^7/nV = 0- (23) 

In terms of Hirota's ^^-derivatives this can be written as (4), with the function 

P(X, Y):=gs e^+^ + ge e^"^ + + h + h. 

However, since the D-operator is antisymmetric and P operates on a symmetric target 
/ • /, only the symmetric part of P is relevant. Thus the proper form is 

P(X, Y) := ^8 cosh(X + Y)+g6 cosh(X -Y)+g7 cosh(X) + he cosh(y) + hj. (24) 

Since equation (9) still holds, the construction of ISS and 2SS goes as usual, resulting 
with formulae (5-8), with rjj — Pjn + qjm+rjj. In particular, from the OSS fn,m — 1; Vn, m 
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we get the condition 

g^ + gi + gs + h^ + h-!^ o, (25) 

and from the ISS fn,m — 1 + e'^^ we get the DR P{pj, qj) — 0, i.e., 

Qs cosh(|jj + Qj) + cosh(pj - Qj) + g-j cosh(pj) + cosh(gj) + /ly = 0, (26) 

while the formulae for A^j (8) and for the 3SC (11) hold as written. 

The practical problem is the following: How to implement the dispersion relation 
(26) while keeping the phase factor and the 3SC manageable? It seems that in order to 
implement the DR we must rationalize it somehow, but then we will loose the additive 
form of the phase factor and the 3SC. 

As the first step we replace 

and introduce the function 

P\x, y) := {xy)-\g^ {x^ + 1) + {x^ + y^) + 97 y{x^ + 1) + h x{y^ + 1) + 2h7 xy], 
after which the dispersion relation is still P'{p'j, g^) = 0, but 



A. 



and correspondingly in the BSC. This does not look bad, except that the DR is still hard 
to implement. Next we continue with 

so that 

fn,m=l + C,{^) (1^) . (28) 

Now the DR gets a nice form 

y(P,-, Qj) := (P/ - l)-\Q^j - l)-\oi P/ + 2o3 PjQj + 02 Q] + 04) = (29) 

where 

^ _ ^ , I, ^6 = \{-Ox - 02 + 2O3 - 04), 

_ With inverse g% = jy—oi — 02 — 203 — 04), 

^■^-^^Tf /l6 = i(-Ol + 02 + 04), 

tir — 2(01 + 02 — 04). 
However, the simple additive property in A^ is lost, we find instead 

A - — ' Qz-Qj ) /gpH 

\ Pi+Pj ' Qi+Qj ) 

and the 3SC is still much more complicated. 

In the following we use the parametrization (28), which imphes (29,30). 
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3.2 Results 

The search was performed by computing the 3SS subject to the dispersion relation 

oi Pf + 2o3 PjQj + 02 Q^j + 04 = 0. (31) 

The dispersion relation was implemented using LET-rule of REDUCE [13] on the leading 
monomial (in "gradlex" ordering). In REDUCE such rules are evaluated to the end. 

• For the generic case we may assume that Oi 7^ and scale it to Oi = 1, then we can 
use the rule: for all n let pp(n) "2=-2*o3*pp(n)*qq(n)-o2*qq(n) "2-o4; The 

resulting 3SC is rather long but the coefficients of various powers of Qi and ffist 
powers of Pi are relatively simple and one can proceed by solving them one by one. 

• If Oi = then we may assume 03 7^ and scale it to 03 = —1/2. Then we can solve 
for Pj and use the substitution for all n let pp(n) = (o2*qq(n) ~2+o4)/qq(n) ; 
Some of these cases are n m reflections of Oi = 1, 02 = 0. 

The process branches somewhat but when the results are collected and sub-cases 
eliminated we get the following cases: 

1. oi = 1, 02 = 203 — 04 - 1 

/n+l,m+l/n-l,m-l03+/n+l,m/n-l,m(03— 04 — l)+/n,m+l/n,m-l( — 03+l)+/n,m(~03+04) = 0. 

2. Oi = 1, 02 = — 2o3 — 04 — 1 

/n+l,m/n-l,m(03 + 04 + 1) + fn+l,m-lfn-l,m+l03 - fn,m+lfn,m-l{o3 + 1) - /n,m(o3 + O4) = 0. 

3. oi = 1, 02 = 1, 04 = 

fn+l,m+lfn-l,m-l{03 + 1) + /n+l,m-l/n-l,m+l (" O3 + 1) " '^fn,m = 0- 

4. oi = 1,02 = 0,04 = 1 

/n+l,m+l/n-l,m-l(03 + 1) " '^fn+l,mfn-l,m + /n+l,m-l/n-l,m+l(— O3 + 1) = 0. 

5. 01 = 0, 03 = -1/2, 04 = 02 

/n+l,m+l/n-l,m-l( — 2o2 + 1) + /n+l,m-l/n-l,m+l (" 2o2 — 1) + '^fn,m+lfn,m-l02 = 0. 

6. oi = 0, 03 = -1/2, 04 = 1 - 02 

/n+l,m/n-l,m(— 2o2 + 1) — fn+l,m-lfn-l,m+l + fn,m+lfn,m-l + /n,m(2o2 — 1) = 0. 

7. Oi = 0, 03 = -1/2, 04 = -1 - 02 

fn+l,m+lfn—l,m—l ~^ fn+l,mfn—l,mi 202 1) fn,m+lfn,m—l 

+ /n,m(202 + l) = 0. 

For further analysis we restore the freedom in the scaled parameters and combine cases 
to the extent possible. The final results can then be divided into the following cases: 
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Case lb) Case 2) Case 3a) 

Figure 2: Points involved in the equations with 3SS. Case la is obtained by a 90° rotation 
from Case lb, and Case 3b the same way from Case 3a. 

Case 1: 01 + 02 + 04 = ±203: 

fn+l,m+lfn-l,m-l03 + fn+l,mfn-l,m{02 — 03) + fn,m+lfn,m-li0l ~ Os) + fn,mi^3 ~ Oi — O2) = 0. 

(32a) 

fn+l,m—lfn—l, m+lOs — fn+l,mfn-l,m{02 + O3) — fn,m+lfn,m-l {Ol + O3) + ^(03 + Oi + O2) — 0. 

(32b) 

These are related by m ■<-> — m; due to (28) this corresponds to Q — Q, or in the 
dispersion relation, 03 ^ — 03. The associated diagram is given in Figure 2. 

For Ol = 1 this contains items 1 and 2 of the previous list, and for Oi = 0, 03 = 
— 1/2, 04 = — 1 — 02 the first one also contains item 7 and for oi = 0, 03 = — 1/2, 04 = 1 — 02 
the second one contains item 6. 

Case 2: oi = 02 7^ 0, 04 = (item 3 from the previous list): 

/n+l,m+l/n-l,m-l(c3 + Oi) + /„+l,m-l/n-l,m+l (" O3 + Oi) — 20i/^^^ — 0. (33) 

This is the Toda-lattice (center + cross, see Figure 2). 

Ccise 3: 02 = 0, 04 = Oi or oi = 0, 04 = 02 (or: 01 + 02 = 04, Oi — 02 = ±04). 

fn+l,m+lfn-l,m-l{03 + Oi) + fn+l,m-lfn-l,m+l{oi — O3) — 20ifn+l,mfn-l,m = 0. (34a) 
fn+1 ,m+l fn—l,m—l (02 + O3) + fn+l,m-lfn-l,m+l{02 " O3) - 2o2fn,m+lfn,m-l = 0. (34b) 

These contain item 4 and its n ^ m reflection, item 5. Both are three term equations with 
two diagonals and a vertical or a horizontal, see Figure 2. They correspond to Hirota's 
discretization of the KdV equation [16]. 

4 Post processing 

In this section we will take a closer look on the equations obtained, and in particular on 
how the soliton solutions look and how the final results can be obtained as a reduction 
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from the Hirota or Miwa equations. In the following the parameters o, are from the 
search, a, b, c are those used in the Hirota and Miwa equations (17a), (18a), and a, (3 are 
the canonical parameters in terms of which the PWF's have the simplest form. 



4.1 Case 2: 

4.1.1 A better pcirameterization 

The dispersion relation for Case 2 reads 

oiP^ + 2osPQ + oiQ^ = 0. (35) 

Parameterizing 01(7^ 0), 03 as 

Oi^2ap, 03 = -(q;2 + ^'), (36) 

the DR factorizes as 

(aP - I3Q) {I3P - aQ) = 0. (37) 

Let us fix a, (3 as determined from (36). In the generic case 7^ (3'^ and we have in fact 
two different kinds of solitons, depending on which factor of (37) is satisfied. Thus we 
may choose either 

Qj = ^/kj, Pj = a/kj, 



or 



and obtain the PWF 



Qj = a/kj, Pj = /3/kj, 



The phase factor Aij has the canonical form 

^"-(Mf- '''' 

4.1.2 Reduction from Hirota's equation 

Equation (33) is obtained from (17a) by the reduction 

T(n+ l,m+ 1,A;) = r(n, m,A; + 1), (40) 
which we use to define a new function / by: 

Thus we define the new independent variables by 

n' :— m -\- n -\- 2k — 2, m' :— n — m 
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and then it is easy to see that the 6 corner points of the cube on which (17a) is defined, 

namely {(100), (Oil), (010), (101), (001), (110)} will project to {(-1, 1), (1, -1), (-1, -1), (1,1), 
(0,0), (0,0)}, displayed at the center of Figure 2. With this definition one obtains (33) 
for fn'm' with, 

oi = |c(6 — a), 03 = |(c(a + b) — 2ab). 



4.1.3 Reduction of solution 

For the NSS of DAGTE given in (17b) the reduction (40) imphes 



1 — aqi\ f 1 — hq, 



_ 1 - cqi 
I - apij \l-bpi) \l- cpi 

which has the nontrivial solution 

ahpi — a — b + c 



ab{cpi - 1) 



(41) 



Here we must assume that ab{a — c)(6 — c) 7^ 0. After this the PWF of (17b) gets the 
form 

1 N n'+2 / , , , 1 ■ 

1 \ f rrz r^bpi-1 



' a—c I b—c 
b 



" CPi-1 

We can write this in the canonical form for 2D systems 

\a + kj \p + ki 
where instead of p we have used k defined by 
1 — uv ki + s{u — v){l + uv 



a —c 
b V b-c 



api - 1 



(42) 



(43) 



Pi 



c ki + s{u — v){l — uv) 



where u :- 



b ' 



V :- 



(44) 



and then 



a — s{uv — l){u — v), (3 — s{uv + l){u + v). 



Here s is an arbitrary parameter. Note that for this case we have already derived PWF 
in the form (38), which allows different kinds of solitons by a •<->■ /3. In the above the sign 
change in u corresponds to this freedom. 



4.1.4 Special cases 

The case Oi = 03 corresponds to a{b — c) = and Oi = — 03 to b{a — c) = 0, both of which 
break the formulae in Sec. 4.1.3. In fact in these cases the equation (33) only has two 
terms and the system can be reduced to the ID equation f^+ifN-i ~ In = ^■ 

The case 03 = corresponds to a = ±i^, which imphes c = 2ab/ (a + 6). The above 
formulae still hold for this special case. 



13 



4.2 Case 3a: 

4.2.1 A better parameterization 

For Case 3 the DR is (01O3 7^ 0) 

OiP-^ + 2o3Pigi + oi = 0, i.e., Q,^--^{Pi + l/Pi), 

203 

after which the PWF of (28) becomes 

(P.-lY ( o,Pf + 2o,P, + o, y ^ _ / (P.P,-1)(P.-P,^ ^^ 

\P^ + l) \0iP,^-20sPi + 0j ' \{P,Pi + l){Pi + Pj 

In order to simphfy this we define a new sohton parameter ki by 

2aP 



ki 



and furthermore let 03 = —aoi/f3, after which we obtain 
with Aij as in (39). 

The appearance of n/2 as the exponent in (45) means that the natural equation in 
this case is obtained by replacing n — 2N, i.e., 

/Ar+l,m+l/iV,m-l(03 + Oi) + /Ar+l,m-l/Ar,m+l(oi — O3) — 20i/iv+l,m/Ar,m — 0, (46) 

which involves points in the 2x3 sub-lattice. (In fact this corresponds to the three-term 
equation (20) in [20].) 

4.2.2 Reduction of equation 

Equation (34a) is obtained from Hirota's equation (17a) by the reduction (see also [17], 
Part III, Example 1) 

T{n,m -\- l,k + 1) — T{n,m,k), (47) 
which we use to define a new function / by: 

'^n,m,k f2n,m—k- 

Defining further the new independent variables by 

n' := 2n — 1, m' := m — k, 

one obtains (34a) for / with, 

Oi = |a(c — 6), 03 = l{a{b + c) — 26c). 

In this case the 6 corner points used in (17a) are projected to {(1, 0), (—1, 0), (—1, 1), (1, — 1), 
(—1,-1), (1,1)}, as displayed in rightmost part of Figure 2. 
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4.2.3 Reduction of the solution 

As for the PWF in (17b) the reduction (47) means 



^-bpij \l- cp, 



= 1 



which has the nontrivial solution 



1 1 

qi = -Pi + T + -■ 

c 

After this the PWF of (17b) becomes 

'a{b + c) - bc{api + i)\ ("'+i)/2 / _ cp^ 



bc{api - 1) 7 vc(-l + bpi) 

We can write this in the form (45) after redefining p ^ khy 

where we have used the parameterization 

/ n ■,\ 1 a + (5 + 1 
a — cya -\- p -\- 1) ^ b — c — -, oi — —pbc, o^ — abc. 

a — p + 1 

The phase factor Aij is as in (39). 
4.2.4 Special cases 

The above construction works except if Oi(oi + o^){oi — 03) = 0, corresponding to abc{a — 
b){a — c){b — c) = 0. In these cases equation (34a) reduces to a simple 2D equation, which 
can be normalized, e.g., to /jv,m/jv+i,m+i - /w+i,m/w,m+i = 0. 

4.3 Case lb: 

We consider only Case lb, i.e., 04 = — (oi + 02 + 203), Case la is similar with some sign 
changes. Here we may assume that 03(03 + Oi)(o3 + 02) 7^ because otherwise the system 
reduces to one of the 3-term equations discussed already. 

4.3.1 Details 

The dispersion relation is given by 

oi + 203 PjQj + 02 Q] = 01 + 02 + 203. (48) 
In the generic case this can be resolved by 

Pi = TT^(% + Q3 = TT^(% + ^/%), (49) 

15 



with 



Oi = 2-i(3{l + a)\ 02 = 27a(l + /3)2, 03 = -7(a + + a)(l + = 2-i{a- (3f . 

(50) 

Here we must assume that + 7^ 0. Furthemore, the three special cases {a — (3), 

(a = 0, /3 = 1) and (a = 1, /3 = 0) may also be omitted as they lead to two-term equations. 
The PWF now becomes 

^ Uki-l){ki-a) 
[{ki + l){ki + a) 

and the phase factor 

klklja + f3-l)-{kl + ki)aP + kik2{a + I3){a- - 1) + a(3{-a/3 + a + /3) 
klkl{a + /3 - 1) - (A;? + kl)a^ - kik2{a + /3){a- l){/3 - 1) + a/3{-a/3 + a + /3)' 

Note that a,/3,7 as determined from (50) are not unique, but rather one can change 
ai— )■ 7> 71-^ ja'^f]'^. If this is accompanied with ki i->- l/Zcj then equations 

(51a,51b) do not change. 



(h + l)(h + /3) 



(51a) 



4.3.2 Reduction from Miwa's equation 

For Case 1 the equation has four terms and therefore one should be able to obtain (32b) 
as a reduction of Miwa's equation (18a). Let us consider the reduction 

r(n, m, k) = T{n + 1, m + 1, A; + 1), (52a) 

and define 

'^n,m,k — fn—k,k—m (52b) 

and change variables by 

n' = n — k, m' = k — m, (52c) 

This projection takes the corners of the cube on which (18a) is defined to the points given 
on the left of Figure 2, in particular the diagonal {(0, 0, 0), (1, 1, 1)} projects to (0,0). 
Then for the equation we obtain (32b) for n', m' (we drop the primes in the following) 
with the parameter correspondence 

oi = 2a{b^ - c^), 02 = 2b{c^ - a^), 03 = (a - b){a + c)(6 + c), 04 = -2c(a^ - 6^). (53) 

Here we may assume that (a^ — 6^) (6^ — c^)(c^ — a^) ^ 0, because otherwise we get a 
two-term equation (22) with Ci + C2 = 0. 

The solution to Miwa's equation (18a) is given through (18b, 18c) and the reduction 
(52a) implies the condition 

(1 - api){l - aqj) (1 - bpi){l - hqj) (1 - cpi){l - cqj) _ ^ 
(1 + api){l + ag^) (1 + hpi){l + hqi) (1 + cpi){l + cqi) ~ 
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and since Pi + Qi this leads to the bi-quadratic 

abc{ab + ac + bcjp'^qf + abc{p1 + q^) + (a + + c)(c + a)piqi + a + 6 + c = 0, (54) 

which suggests an eUiptic parameterization (if ahc ^ 0). This form of solution agrees with 
(51a), if in addition to (54) we have 

{ki - l){ki -a) api){l - aqi) ^^^^^ 



(/Cj + l)(fci + a) {I + api){l + aqi] 



(55b) 



(fc, + l)(fc, + /3) {l-hp,){l-hq. 

From these 3 equations one should, in principle, be able to determine p,, qi in terms ki 
and a, ^ in terms of a, 6, c, but this does not lead to simple expressions for them. In any 
case we have verified that the equations are consistent provided that a, /3 are determined 
from 

I in 1 ^ b(a — c)'^ ^ a(b — c)'^ . , 

ra + s^^l, r/Q; + s//3 = l, where r := — — r, s := — -^r. (56) 

c(a — by c{b — ay 

The same equations are obtained by equating (50) and (53). Furthermore, one can show 
that the forms for Aij given in (18c) and (51b) agree, subject to relations (55) and (56). 

4.3.3 Special cases 

The discussion above holds in the generic case, where in particular the curve defined by 
the dispersion relation is irreducible. In fact the Jacobian of the transformation (50) is 
given by 

- -47 {a - P) (1 + a) (1 + /3) = 4(oi02 - O3 



a(a,/3,7) 
Similarly for (53) we get 

9(01,02,03) 



a' - b'){¥ - c'){c' - a') = 3oi0204/(a6c). 



9(a, 6, c) 

We need to consider separately the cases where the RHS vanishes. It is perhaps not 
so evident that in the second case the Jacobian also vanishes when 03 + 04 = 0, since 
03 + 04 = (a — 6) (a — c){b — c). Among the various special limits we are only interested in 
those not included in Cases 2 or 3. 
We now consider these special cases. 

• Special case 01O2 = O3, Oj 7^ 0: This is a singular point for the parametrization (50) 
in the analysis of Sec. 4.3.1. 

We use the relation to eliminate 02. Then equation (32b) becomes 

/n-l,m+l/n+l,m-l O1O3 — /n-l,m/n+l,m 03(01 + O3) — /n,m-l/n,m+l Oi(oi + O3) 

+ flm {oi + 0103 + ol) = 0, (57) 
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and the dispersion relation (48) 

{OiPj + OsQjf = (Oi + 03)^ 

which is solved by 

Pj = -o'sQj + aj{l + o'g), o-J = 1, o^j = 03/01. 

Clearly there can be two kinds of solitons, depending on the chosen sign aj. 
By directly solving (57) one finds the PWF 



-o'sQj + aj(l + o'^) + lJ \Qj + l 



(58a) 



and the phase factor 



~ QI - {\a,a^ - DQiQj + Q] - {Q, + Q^){a, + a^){\ + ^) - §(1 + 0,0^) ' ^ 
The results (58) can be obtained from (51) with the singular limit £ ^ using 

OL = -1 + 2£0i + 0(£^), 
^ = -l-2£03 + 0(£'), 
7 = -l/(8£%)+0(£-^), 

and 

hi^Gi- £Oz{Qi - (Ti) + 0(£^), 

where £ is defined by O1O2 — ol — —£^03(01 + 03)^. 

As for the the Miwa's parametrization (53) the relation 01O2 = o| is a regular point. 
For convenience let us use the parametrization 

a = v{l + s)/{?,s-l), 6 = w(l -s)/(3s + 1), c=-sv, 

we get 

(352 + l)(s+l)V _ (3s2 + l)(s-l)V _ (3s2 + l)(s2_l)^;3 

^^-"^ (3^Tl? ' ""'-'^ (3^^1? ' ""'-^ 9i^^l ' 

where s are the new parameters. Then the reduction condition (54) is satisfied by 
choosing 

Gi Mi — 1 (jj 3^1 + 1 

V + 1 V li — 1 

which again shows the possibility of two different kinds of solitons, due to (Tj = ±1. 
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With these choices the phase factor and the PWFs of (18) only depend on £^ and it 
is then convenient to replace 



3A;'-4r' 



2(35^ + l)^(s + iy 
(3s + l)2(3s-l) ' 



(59) 



M 

^k'i + (oi + 203) + SajOi J \k'-- (2oi + 03) - SajOs 
We note that r solves 



and with this we get form (18) the symmetric formulae (recalling (52c) and dropping the 
primes) 

' k'i + (oi + 203) - SajOi \^ (k[- (2oi + 03) + 3(Tj03 ^ ^^^^^ 

(60b) 
(61) 

which appears later as well, and that the transformation (59) is singular if r = 0. Finally 
the result (60) is related to (58) via 

k'- = -3crj03(5j + 2oi + 03. 



k[^ + k[k'2 + k'^^ - 12(of + 01O3 + ol) 



ol + O1O3 + O3, 



• Special case 03 + 04 = 0: (i.e., oi + 02 + 03 = 0). This is a regular point for the 
derivation of Sec. 4.3.1 but a singular point for Sec. 4.3.2. It leads to a 3-term equation. 



O3 /n+l,m-l/n-l,m+l + Ol fn+l,mfn-l,m ~ (Ol + O3) fn,m+lfn,m-l — 0, 

which is different from the 3-term equations in Cases 2 and 3. 

In the analysis of Sec. 4.3.1 we should just parameterize a, ^, 7 by 



(62) 



a = (oi+r)V(o3(oi+03)), /3 = -{r+o,+OsY /{oiOs), 7 = -{2r^-ol-r{2oi+Os)y /{8osr^), 

where is as before (61). If we now reparametrize kj = K.j{2oi + 03 + 2r)/(3o3) then the 
PWFs (51a) become ^ 



X 



{o^Kj + 2oi + 03 - 2r){o2Kj + 2oi + 02 + 2r) 
(o3fi;j - 2oi - 03 + 2r){o2Kj - 2oi - 02 - 2r) _ 
{o^Kj + 2oi + 03 - 2r){oiKj + 203 + Ol + 2r) 
(o3fi;j - 2oi - 03 + 2r)(oifi;j - 203 - Oi - 2r) 



and 74jj (51b) factors as 



A. 



[(^KiKj 3)^ (Kj + Kj)^] (Kj '^j)^ 
[(KjKj- + 3)2 - {Ki - Kj^liKi + KjY ' 



(63a) 



(63b) 



^For these computations it is convenient to use the rational parametrization 0\ = z{^w'^ — l)/{2>w^ + 
1), 02 = [4(3w-l)wz]/[(3w^ + l)(w-l)], 03 = -[(3w-l)2(w+l)z]/[(3w2 + l)(w-l)] and correspondingly 
a = {w + l)/[(3u; - l)w;], /3 = [(3u> + l)(w; + l)]/[(3w - !)(«; - 1)], 7 = z[(3u; - lf{w - l)w^]/[2{iw'^ + 
Vf{w + 1)], r = —{Zw — l)z/{w — 1), where w, z are the two new free parameters. 
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For the singular hmit that should be used with Sec. 4.3.2 we start with 

Ol = -(02 + 03)(1 + 02036^), 



and get the expansions' 



3 



a = l + §(03 - 02)e^ + 0(e6), 
b = l-§(02 + 203)e' + 0(e^), 
c = l + mo2 + Os)e^ + 0{e^), 



and 



2rK{2oi + 03)[3(«;2 - 1)03 - Ai^ar] + 3oscrr{K'' + - 1)03 + Ai^ra] 3 , 6^ 
^^^^ = -'"-^ 3[16«V2 + 3oi(«^-l)^] ^ 

Then in the limit e — > one obtains (63) from (18). 

• Special case oi = u, 02 — oJ^-, 03 = 1, (^^ + (^ + 1 = The combination of the previous 
two cases implies that Oj are cubic roots of unity. For convenience we choose Oj — uj^, 
where cu^ — l,uj ^ 1, and then the equation becomes 

fn—l,m+lfn+l,m—l ~l~ ^ fn—l,mfn+l,m ~l~ ^ fn,m—lfn,m+l 0- 

For analyzing this case we have to take a singular limit in both methods. We approach 
the limit with oi = a; + e^, e — >■ 0, and then we get the expansions {e — iou^e) 

a = -l-2iuje^ + 2uj'^e^ + 0(e^), 

P = -l + 2ie^ + 2e^ + 0{e^), 

7 = luh-'-il(u + 2)e-' + 0(l) 

and (directly or as limits of the 03 + 04 = case) 

a = l + |(a; + 2)e3- |(3a; + 4)e^ + 0(e^), 

b = l + |(a;-l)e3-|a;e^ + 0(e9), 

c = l-|(2cu + l)e^-i(a; + 4)e^ + 0(e^). 

Next in order to have nontriial PWFs and relations (55) we must also take singular hmits 
in k,p,q, namely 

kj = aj[l + qe'{-^ + 2{u;-l))], 

P3 = a' + 2'^{a,{l-uj) + a'{co + l)) + 0{e'), 

qj = -a' + 2'^{a,{l-u)-a'{u+l)) + 0{e'). 



^This is also the type of limit by which Miwa's equation reduces to the Hirota's equations. 
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Here aj,a' — ±1. In the limit e — > we then get from (18) 

''''' " \ky-2uj-l + 3aj ) \k'j - 2uj - 1 - 3aj J ' ^^^"^^ 

and 

(k' - k'V 

^ p-^ ^ il , (64b) 

which also follow directly form (60). These also follow as singular limits from (63), if we 
use Kj = —aj{2u + 1)(1 + ie^kj/S), along with oi — u + e^, 02 — uj^ — |(2a; + 3)e^, 03 = 

-Oi-02 + 0{€^). 

• Special case Oi = 0, O2O4 ^ 0: This can be obtained from the analysis of 4.3.1 by 
just putting ^ = but for (53) a singular limit is needed. Thus we let 

a = 04/02, 13 = 6^04/(403), 7 = 02/(204), 

a = am/e, b = -e^ 02/(2a^), c = -e^ 04/(2a^), where = -203(02 + 03). 

Then with the parametrization 

_ -4(o2 + 03)03 _ -(04 + 02fc,^)e 

Q-m02^'^ ki 2amOzki 

we find that in the limit e ^ equation (53) is satisfied and the formulae of Sec. 4.3.2 
work. Prom both approaches the result is 

^^_ nk,-i){K-a) Y (K-iy _ {k, - k2f (fcifc2 - ^) , . 

\{h + l){h + a)) \h + l) ' {k^ + k2f{k^k2 + ay ^ ^ 

The difference equation in this sub-case reads 

(/n+l,TO-l/n-l,m+l ~ fn,m+lfn,m,-l)o^ — {fn+l,mfn-l,m ~ fn,m)i^3 + O2) = 0. (66) 

The case Oi = 0, O2O4 7^ is similar. 

The further specialization to Oi = 04 = can be obtained from the above by letting 
first 02 = — 203 + pe, i.e., 04 = —pe — e^. This corresponds to a ^ in (65) and 

fn+l,in—lfn—l,m+l fn,m+lfn,m-'l fn+l,mfn—l,m fn,m ^' l^"^) 

Similarly the case Oi = 02 = can be obtained by further putting 02 — e. This 
corresponds to a ^ cxo, which yields 

and 

fn+l,m—lfn—l,m+l fn,m+lfn,m—l fn+l;mfn—l;rn ~l~ fnm ^' l^'^) 

Note that here A12 is negative so this system can have a non-singular 2SS only if it 
composed from a complex conjugate pair of parameters k^. In both of these double zero 
cases the solitons have a fixed angle in the n, m-space. 
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• Special case 04 — (oi + 02 + 203 = 0): This is a regular point for Sec. 4.3.2 but a 
singular point for Sec. 4.3.1. The dispersion relation is 

{P~Q){o,P-O2Q) = 0, 

and thus wc have two kinds of solitons if 02 7^ Oi. The nontrivial factor {oiP — 02Q) is 
included in the parametrization (49) provided that /cj — ?■ together with 04 — )■ 0. We will 
therefore take e := 04, ki = k[e/{oi — 02Y with k[ arbitrary (the (oi — 02)^ factor is for 
simplicity later ) . For the remaining parameters we construct a power series expansion in 
e: 

\0\ — O2) [Ol — O2) 

The limit of the PWF is given by 

the phase factor gets the form 

''~\K + k'J \020,-k[k'J' ^''^ 
and the bilinear equation (32b) reads 

(02 +Ol)/n+l,m-l/n-l,m+l + (^2 " Oi) fn+l,mfn-l,m ~ (^2 — Oi) fn,m+lfn,m-l ~ (^2 +0i)/^ m ~ ^- 

(73) 

In order to keep the connection to the Miwa equation as given in Equations (53-56) 
we need first of all the expansions for a, b, c. The leading terms are 

, Ol Ol 1 s On 

a = ao + . . . , = — ao h . . . , c = qq-^ + . . . , where an = - — . 

02 of — O2 2oi 

Finally to connect p,, qi in (54-55b) to k'^ we find for pj, qi the expansions 

_ 02 _ 2 (oi + 02)(oi - 02falkl 

aoki e {k[ -0102)0^ 

It is easy to see that this satisfies (55) but in order to satisfy (54) is actually necessary to 
compute the second order in the e-expansions, while for (56) three orders are necessary. 

• Special case 04 = with others. The special case where 04 = together with 
02 = Ol, or Ol = 0, or 02 = leads to dispersion relations Q{P — Q) = and P{P — Q) = 0, 
respectively, which are perhaps too simple to be of interest. 



22 



5 Conclusions 



Hirota's direct method has turned out to be very efficient in deriving sohton solutions for 
a given equation. It has also be used as a method for searching for integrable equations 
in the continuous case, the criterion being the existence of three-soliton solutions. 

Here we have applied the three-soliton condition to two-dimensional partial difference 
equations. In comparison with the continuous case there is now much more freedom, since 
the same partial derivative can be discretized in many ways. In this paper we restricted 
the search to partial-difference equations defined on a 3 x 3 stencil. For the 2 x N case, 
see [20]. 

The master bilinear PAEs arc Hirota's three-term DAGTE equation (17a) and Miwa's 
four term equation (18a). In the detailed analysis of the results found here we showed 
that all of them can be obtained from either one of these master equations with a suitable 
projection. However, sometimes a nontrivial singular parameter limit was also needed. 
Since the phase factor of the soliton solution is the same for a discrete equation and its 
continuum limit we can say that most of the results can be considered as discrete versions 
of the KdV equations, but some are related to the Boussinesq or Ito equations, while some 
cannot yet be definitely identified. 
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